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Abstract 

In this work, we have used the hquid-drop model in the context of stabihzed 
jelhum model, to study the stability of Z-ply charged metal clusters of different 
species against fragmentation. We have shown that on the one hand, singly ion- 
ized clusters are stable against any spontaneous fragmentation, and on the other 
hand, the most favored decay process for them is atomic evaporation. However, 
multiply charged clusters of sufficiently small sizes may undergo spontaneous decay 
via fission processes. Comparing the results for different species show that for fixed 
N, the lower electron density metal clusters can accommodate more excess charges 
before their Coulomb explosions. This comparison also shows that, for fixed Z, the 
atomic evaporation which is the most favored decay mechanism for sufficiently large 
clusters, takes place at lower A^s for lower electron density clusters. 
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1 Introduction 



Study of the stability of charged metal clusters against fragmentation is important in 
nano-sized systems. In this work, we have studied the binary fragmentation of A^-atom 
Z-ply charged A4ff clusters in the framework of the liquid-drop P] model (LDM) with 
stabihzed |2I jeUium model (SJM) coefficients for the energies. Here, Z=l, 2, 3, 4; and 
Ai covers the species Al, Ga, Li, Na, K, and Cs. Here, in the case of Al, when A^ is a 
multiple of 3, it corresponds to real Al^ cluster with u = N/3. The clusters are assumed 
to be spherical with sizes determined hj R = N^/^rf with A^ < 100. The quantity rf is 
the Wigner-Seitz radius of the electrons in the bulk metal. The possible decay channels 
oi M% are 

M% ^ Mlzf + M'^^ (1) 

Here, p is positive integer as A^, and < Zi < [Z/2]. Different processes are defined by 
different values of Zi. The process in which Zi = (one of the fragments is neutral), is 
called evaporation. The processes in which < Zi = Z/2 are called charge-symmetric 
fission, and the charge-asymmetric fission is defined by < Zi < [Z/2]. For any given set 
of values A^, Z, Zi, different values of p define different decay channels. In a particular 
decay channel of the evaporation processes, the negativity of the difference between total 
energies before and after the fragmentation, 

D^{N,p) = E^{N -p) + E^{p) - E^{N), (2) 

is sufficient for the decaying in that particular channel. In the above equation, E^{N) and 
E^{N) are the total energies of Z-p\j ionized and neutral A^-atom clusters, respectively. 
However, in fission processes a negative value for the difference energy is not a sufficient 
condition for the decay of the parent cluster. This is because, the competition between 
the short-range surface tension and the long-range repulsive Coulomb force may give rise 
to a fission barrier. The heights of the fission barriers are calculated using the two-spheres 
approximation In Fig. Q the fission of a Z-plj charged A^-atom cluster into two clusters 
of respective sizes Ni, N2 = N—Ni and respective charges Zi, Z2 = Z—Z\ is schematically 
shown. Q/ is the energy release, B^. is the fusion barrier (or Coulomb barrier) which is 
the maximum energy of the Coulomb interaction of two positively-charged conducting 
spheres, taking their polarizabilities into account. -B/ is the fission barrier height which 
is defined as 

Bf = -Qf + B,. (3) 
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The Coulomb interaction energy, Ec, as a function of their separations, d, for two 
charged metaUic spheres can be numerically calculated using the classical method of 
image charges [S]. 

The most favored decay channel in evaporation processes is defined as the channel for 
which the dissociation energy assumes its minimum value 

D^{N,p*) = mm [d^{N,p)} , (4) 

and the most favored decay channel in fission processes is defined as the channel for which 
the fission-barrier height assumes its minimum value, 

Bf{N,p*)=mm{Bf{N,p)}. (5) 
The energy of an A^-atom Z-plj charged cluster in the LDM is given |31 13] by 

E-(7V) = E»W + Z(H- + £) + ^. (6) 

in which W, c, R, and e are the work function of the bulk metal, the finite-size correction 
to the work function, the radius of the cluster, and the electron charge, respectively. For 
simplicity, the position of the centroid of excess charge, a, is neglected in our calculations. 
E^{N) is the energy of a neutral A^-atom cluster in the LDM, which is given by 

E\N) =eN + An{rffaN^/^ + 2Tir^-fN^'\ (7) 

In Eq. ((Zj), the quantities e, a, and 7 are total energy per electron of the bulk, surface 
energy, and curvature energy, respectively. We have calculated the quantities c, W , a, 
and 7 by fitting to the self-consistent Kohn-Sham jE] results in the SJM for different 
values However, e is calculated using the SJM energy expression for the bulk system 
(Eq. (1) of Ref. [7j). 



2 Results and discussion 

Using the method of image charges, we have calculated the Coulomb interaction energy of 
two charged metallic spheres, taking their polarizabilities into account. The calculations 
show that the maximum of the interaction energy, is achieved for a separation > 
Ri + R2- Fig. |21^a) shows the Coulomb interaction energy of an A^i-atom cluster with 
another A^2-atom of respective excess charges Zi = 3 and ^2=1. The radii of the clusters 
are calculated from R = N^^^rf, as in the SJM. In this figure, we have taken r^=2.07 
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which correspond to bulk Al. When both A^^i and N2 are multiples of 3, the results 
correspond to real Al clusters. The value at the maximum specifies the quantity -Be, and 
the position of the maximum, do, is the separation between the centers of the two spheres. 
As is seen, for fixed values of charges. Be is the highest when the sizes are equal. In Fig. 
I2fb), the situation is shown for equal charges Z\ = Z2 = 2. In this case, when the sizes 
are equal, is maximum as before but here, do = Ri + R2. An other feature shown in 
Fig. |2^c) is that, when both charges and sizes are asymmetric. Be is higher if the smaller 
charge corresponds to the smaller cluster. In Fig. I2td), we have compared the Coulomb 
interaction energies of two Al metallic spheres with respective sizes Ni = 2 and N2 = 18 
for different charges. It is seen that keeping the sizes fixed but increasing the charge on 
any one of them, without changing the charge on the other one, increases the height of 
the maximum. 

In our calculations we need the values of the Coulomb interaction energy at the maxima 

points, i.e., the BcS. In Fig. Efa), we have compared the Be values for equally charged 

1/3 

Z\ = Z2 = 1, Al metallic particles but with different radii R\ = 2.07A^i and R2 = 

1/3 

2.07A^2 . This fi gure shows that when both of the sizes are small, the barrier is higher 
than the case when at least one of them is larger. In Fig. E^b), we have compared the 
-BcS for singly ionized pairs Al, Na, and Cs with different values of A^^i and A^2- In the 
aluminum case, since the radii of the spheres are smaller than those of Na and so on, their 
polarizabilities are thereby smaller and thus the Al values lie above that of Na, and so on. 
Now, if the charge on one of the clusters is increased without changing that of the other 
one, the arrangement of the curves would not change but they shift upward according to 
Fig. Hd). 

In the energy calculations of Eq. © we use the results obtained for c and W in Ref. 
IS). The values of e are calculated using Eq. (1) of Ref. [7j for unpolarized case, a 
and 7 in Eq. (|7j) are obtained by fitting the results of the self-consistent solutions of the 
Kohn-Sham equations in the SJM for neutral clusters of sizes < 100, and different 
values. The values of rf for Al, Ga, Li, Na, K, and Cs are 2.07, 2.19, 3.28, 3.99, 4.96, 
and 5.63, respectively. In Figs. Ill^a)-(c) we have plotted the SJM values of e, a, and 7, 
as functions of r^, respectively. 

To show the variations of the dissociation energy, D^{N,p), for different charging and 
different evaporation channels, we have plotted, in Figs. El^a)-(b), the quantities D^~^{N,p) 
and D'^^{N,p), respectively, as functions of the neutral fragment size, p, for Li. As is seen 
in Fig. Efa), for small p values D^^[N,p) has an increasing behavior, then after passing a 
maximum, it changes to a decreasing behavior, and finally, for large p values it increases 
again. This behavior has its roots in the trade off between the surface energy (the volume 
energy term does not change in the fragmentation process) and the last term of Eq. @. 
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For small p, the surface difference 



A5(i?i, R) = Att [rI + {R^ - Rlf'"^ - R^ 

increases, and then for intermediate p values, it assumes a maximum, and finally decreases. 
The increasing behavior at the tail is because, with increasing p, the size of the singly 
ionized fragment decreases and hence, beyond a certain p value, the last term in Eq. (jU)) 
dominates the surface energy term. However, beyond a Z value for the parent cluster, 
the last term of the Eq. © dominates the surface energy term for smaller p values before 
D^{N,p) reaches its maximum as in the case of Fig. El^b). In fact, we have found that for 
Li^Q this transition takes place at a point 1.6 < Zt < 1.7. The value of Zt depends on N. 
This dependence is clearly shown in Fig. Efa) by the fact that D^~^{20,p) and D^^{10,p) 
are increasing functions of p in the whole range. 

The dissociation energy D^{N,p) for small p/N takes the form 

D^{N, p) 27r7r,p^/^ + 47rar,V/^ (9) 

which specifies the most favored channel in an evaporation process with p* = 1. This 
means that, in evaporation processes, the parent charged cluster prefers to evaporate 
single atoms than larger neutral clusters. This fact is also numerically shown in Figs. 
inta)-(b). On the other hand, the positivity of D^{N,p) in the whole range implies that 
the parent charged cluster is stable against any spontaneous evaporation. 

In Fig. we have plotted the most favored dissociation energies, D^{N,p*), of Na 
clusters as functions of the parent size, N, for different Z values. The results show that 
for large enough clusters, the dissociation energy is independent of Z, which is consistent 
with Eq. ©. However, for small clusters, the dissociation energy increases with charge. 

Figures [7|^a)-(b) compares the most favored dissociation energies of different species 
for Z = 1 and Z = 4, respectively. It is seen that detachment of a single atom from a 
lower electron density cluster is easier, i.e., it needs less energy. Comparing these two 
figures also shows that, although the hierarchy is the same, the dissociation energies are 
shifted to higher energies as we go from Z = 1 to Z = 4. 

Now, we consider the fission processes. We first consider the following two cases: 

Grij^ > Ga^yLp -|- Ga^^, (10) 

Ga^+^Ga^% + GaJ+. (11) 

Processes fllO|) and (fTT|) describe charge-symmetric and charge-asymmetric fissions, 
respectively. In Figs. Efa) and Efb) we have plotted their respective fission barriers. 
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Bf^'{N,p), for Zi = 2, 1 as functions of p, for different sizes of the parent. As is seen in 

z z 

both of these figures, there exists a minimum size for the parent, iV^^„\ beyond which the 
fission barrier height is positive in all possible channels of that particular process. That 

7 Z Z 

is, the parent cluster A4fj with > A'min^ is stable against any spontaneous fission via 
that particular process. 

One of the significant differences in the charge-symmetric and charge-asymmetric fis- 
sion is that [as seen in Figs. Efa) and|Hl^b)], in the most favored channel, the products 
have more or less the same sizes in the former case; whereas, in the latter case the less 
charged fragment assumes quite smaller sizes. That is, the charge-symmetric fission pro- 
ceeds mostly via mass-symmetric channel, and the charge-asymmetric fission does it via 
mass- asymmetric channel. In Fig. Efc), we compare the most favored sizes, p*, for the 
charge-symmetric and charge-asymmetric fission of Na^. It is seen that in the charge- 
symmetric case, p* increases with whereas, in the charge-asymmetric case it remains 
more or less constant. 

To decide whether a cluster with a given size and charge, A^^, is stable against any 
spontaneous decay, one should compare A^ with the quantity 

N^,^ = Tcmx{N^i'}. (12) 

If A^ > A'^in' then the cluster is stable against any spontaneous decay. To show the 
values of Aj^in^ schematically, we plot the most favored fission barriers, B'^'^^[N ,p*) for 
all possible Zi values. The evaporation processes need not be considered because, as was 
shown before, in the LDM there is no spontaneous decay via evaporation. However, in 
order to determine at what size which process is dominant, we include the most favored 
evaporation process as well. It should be mentioned that, in the self-consistent Kohn- 
Sham results [S] , because of the shell effects, the spontaneous evaporation is also possible. 

In Fig. El we have compared the most favored fission barriers, B^'^^{N,p*), of Ga^^ 
for A^ < 150 and Zi = 1,2 with the most favored dissociation energy, D'^^N^p*). We 
have labeled the intersection points of the curves with each other and with the horizontal 
A^ axis by A, B, C, D, and E. At point A, the charge-symmetric and charge-asymmetric 
fission start their competition; thus we name the corresponding size as A^s-a- For clusters 
A^j^ with A^ < Ns_a, the charge-symmetric fission is dominant whereas for larger sizes 
charge-asymmetric one dominates. We name the size corresponding to the point B as 
^min because for sizes greater than this there is no charge-symmetric spontaneous fission. 

4 2 

However, spontaneous fission via charge-asymmetric process still persists for N^^^ < N < 
N^l^. Beyond the size N^]^ which corresponds to the point C, the spontaneous charge- 
asymmetric fission also stops. Our results show that, for small values, A'mrVi — ^mln 
whereas, for > 4.5, the reverse inequality is at work (as in the cases of K and Cs). 
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According to Eq. (fT^ . if > N^^^, then any spontaneous fission stops and the 
cluster A/i^ would be stable against any spontaneous decay. At sizes N^^^ and A'g^^ 
which correspond, respectively, to the points D and E, competition of the evaporation 
with charge-symmetric and charge-asymmetric fission starts, respectively. Beyond the 
size 

A^et = nmx{^f.f }, (13) 

the dominant decay process is evaporation. Our results show that, at least for < 7, 
the inequality A'g^^ < N^^^ is satisfied. When Z = 3, the decay processes consist of 
evaporation and charge-asymmetric fission whereas, for Z = 2, the fission is a charge- 
symmetric one. 

In Figs. ITUT a) and (b), we compare the most favored fission barrier heights B^'^ {N , p*) 
and B'j'^{N,p*) for different species, respectively. The figures indicate that, for large 
enough clusters, the fission barriers are higher for higher electron density metals. 

In Fig. [TTfa), we compare the smallest stable sizes, N^^^ for different species. Looking 
at the figure, one concludes that, for fixed A^, the higher electron density metal clusters 
(say Al) has less capacity for charging than the lower density metal clusters (say Cs). 
This behavior is observed for all 1 < Z < 4. 

To show the size beyond which evaporation dominates all decay processes, we plot 
^eva "with Z = 2,3,4 for all species in Fig. [TTfb). It is seen that, for fixed charge, Z, 
the evaporation process dominates at smaller A^ values for lower electron density metal 
clusters. The trends in Figs. HU^a) and lTTT b) are seen to be the same. 

Finally, in Fig. ITTT c). we compare the Ns_a, the size at which the competition between 
the charge-symmetric and charge-asymmetric fission starts. This size decreases from Al 
to Li, and then increases from Li to K, and finally decreases from K to Cs. In the plot, 
we have taken into account the decimals for the intersection points which can be rounded 
to the next nearest integer. 

3 Summary and conclusion 

In this work, we have studied the stability of Z-plj charged {Z = 1, 2, 3, 4) metal clusters 
of species Al, Ga, Li, Na, K, and Cs using the liquid-drop model with stabilized jellium 
model energies. Fragmentation of clusters into smaller ones can proceed via evaporation 
or fission. Our results show that, in the LDM, any charged cluster is stable against 
spontaneous evaporation, and the most favored channel of induced fragmentation depends 
on the size and the charge of the parent cluster. Sufficiently small multiply charged 



7 



clusters, however, may undergo spontaneous fragmentation via different fission processes. 
We have obtained, for each species, the smallest size that a Z-plj charged cluster is 
stable against any spontaneous fragmentation. We have also shown that, for sufficiently 
large clusters the most favored decay channel is atomic evaporation. Comparing the 
results for different species show that, for fixed N, the lower electron density clusters 
have higher capacity for charging. Results also show that, for fixed amount of charge, 
atomic evaporation dominates at lower A^s for lower electron density metal clusters. 
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Figure 1: Fission barrier in the two-spheres approximation. The parent A^-atom Z-p\y 
charged cluster decays into two clusters of sizes A^i and N — Ni, with charges Zi and Z2, 
respectively. 

Figure 2: Coulomb interaction energy of two charged metal spheres, in electron volls, as 
function of distance. The value of the maximum is Be, and the location of the maximum 
is at do. (a)-the spheres have charges 3 and 1, (b)-the spheres have equal charges Z — 2, 
(c)-the charges of the two spheres are exchanged, and (d)-the sizes are kept fixed but the 
charges are changed. 

Figure 3: Coulomb barrier height. Be, in electron volts, (a)-for two equally charged but 
different sized pairs of the same species, (b)-comparison of Be for different species of 
equally charged pairs with one of the pair sizes kept fixed. 

Figure 4: Functions of r^, (a)-total SJM energy per electron in the bulk metal, (b)-the 
surface energy obtained from fitting to the Kohn-Sham results, and (c)-the curvature 
energy obtained from fitting. 

Figure 5: Dissociation energy, in electron volts, for different evaporation channels and 
different parent cluster sizes. (a)-singly ionized, (b)-4-ply ionized Li clusters. 

Figure 6: Na dissociation energy, in electron volts, as function of the parent cluster size, 
for the most favored channel of different evaporation processes 

Figure 7: Most favored channel dissociation energies, in electron volts, of different species 
for (a)-singly ionized clusters, (b)-4-ply ionized clusters. 

Figure 8: Fission barrier heights of 4-ply ionized Ga clusters, in electron volts, for different 
parent sizes and fission channels of (a)-symmetric, (b)-asymmetric processes. 

Figure 9: Barrier energies, in electron volts, in the most favored channel of different sizes 
of 4-ply ionized Ga clusters. Zi — 1,2 correspond to asymmetric and symmetric fission, 
respectively, while the Zi — case correspond to evaporation process. 

Figure 10: Most favored fission barrier heights, in electron volts, of different species for 
(a)- 3-ply ionized, and (b)-doubly ionized clusters. 
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Figure 11: Functions of Vg, (a)-the minimum stable size for different charges, (b)-the sizes 
at which evaporation dominates fission, and (c)-the sizes at which the symmetric and 
asymmetric fissions start their competitions. 



11 



< ► 

do 




7.5 




8 9 c/ 10 11 12 13 14 15 16 
^ 



d (a.u.) 



12 



11 - 



c\r 9 



8 - 



7 - 



Al, Z1=2, Z2=2 



■0-(N1,N2)=(2,18) 
^(N1,N2)=(5,15) 

O(N1,N2) = (10,10) 




J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



7 



8 



10 



11 



12 



13 



14 



15 



16 



d (a.u.) 




d (a.u.) 



CD 



N 



14 - 



12 - 



10 - 



8 - 



6 - 



4 - 



2 - 



Al, N1=2, N2=18 





0(Z1,Z2)=(1,1) 
^(Z1,Z2) = (2,1) 
♦(Z1,Z2) = (3,1) 
-^(Z1,Z2)=(2,2) 



Q I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

7 8 9 10 11 12 13 14 15 16 17 18 19 20 



d (a.u.) 




N2 




N2 




(a. u.) 



5 




-♦-N=20 



: ^N = 10 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I V ^ 

10 20 30 40 50 60 70 80 

P 



25 




10 20 30 40 50 60 70 80 



P 



CD 




5 10 15 20 25 30 35 40 45 50 55 60 

P 



25 



20 



15 



10 







A/a (rs=3.99) 



^Z=4, Z1=2 
OZ=4, Z1=1 




J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L 

5 10 15 20 25 30 35 40 45 50 55 60 



N 



